clear, close, clc

t = 0:0.02:10;

x = exp(-t) + exp(-2*t);

plot(t,x,'k')
hold on

y0 = [2;-3];
A = [0 1;
    -2 -3];
[t,y] = ode45(@(t,y) A*y, t, y0);
plot(t,y(:,1),'r--')
legend('Analytic','OEDS')
xlabel('Time')
ylabel('Solution x')